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Abstract 

/(ii)-theories of gravity are reviewed in the framework of the matter-antimatter asymmetry in 
the Universe. The asymmetry is generated by the gravitational couphng of heavy (Majorana) 
neutrinos with the Ricci scalar curvature. In order that the mechanism works, a time varying non- 
zero Ricci curvature is necessary. The latter is provided by f{R) cosmology, whose Lagrangian 
density is of the form C{R) ~ f{R)- In particular we study the cases f{R) ~ ii + aR^ and 

f{R) ~ 

PACS numbers: 98.80.-k, 98.80.Cq 
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1. INTRODUCTION 



The observation that the present phase of the expanding Universe is accelerated [Jj] has 
motivated in the last years the developments of many models of gravity which go beyond 
the general relativity, and therefore the standard cosmological model. Among the different 
approaches, the /(i?)-theories of gravity have received a great attention. The reason relies on 
the fact that they allow to explain, via a gravitational dynamics, the observed accelerating 
phase of the Universe, without invoking exotic matter as sources of dark energy. Moreover, 
they also provide an alternative approach Dark Matter problem. 

The Lagrangian density of these models does depend on higher-order curvature invariants 
(such as, for example, R^, R^j^^R^^ , RDR, and so on) j2|, [sj. In this paper we focalize our 
attention to f{R) models which are a generic function of the Ricci scalar curvature R 

^=^2 j d'^^Tg f{R) + S^[g,., ^] . (1.1) 

In particular, the function f{R) we concern are of the forms: f{R) = R + aR^ [n = 2) and 
f{R) = R"", with m = 1 + e ~ C(l) (e < 1). In Eq. ([HI]), is the action of matter and 

= SttG = 8nmp^ [nip c::^ lO^^GeV is the Planck mass). Cosmological and astrophysics 
consequences of (11. ip have been largely studied in literature (see also [2|-0]). 

The aim of this paper is to study the origin of the baryon number asymmetry in the 
framework of f{R) theories of gravity. As well known, the matter-antimatter asymmetry 
in the Universe is still an open problem of the particle physics and cosmology [of. The 



successful prediction of the big bang nucleosynthesis (BBN) [10|, lll[ and the observations 
of CMB anisotropics combined with the large structure of the Universe 

0,1 

I13I I show that 

the baryon to photon number of the universe, i.e. the parameter that characterize such a a 
asymmetry, is of the order 

r^exp = ^^^^ < (9.2 ± 0.5) 10"" , (1.2) 

where Ub {^3) is the baryon (antibaryon) number density, and s = (27r^/45)5'^,r^ is the 
entropy density ((7* ~ 100 are the relativistic degrees of freedom). 

As shown by Sakharov, a (CPT invariant) theory is able to explain the baryon asymmetry 
provided that the following conditions (Sakharov's conditions) are fulfilled ^]: 1) there must 
exist processes that violate the baryon number; 2) the discrete symmetries C and CP must 
be violated; 3) departure from thermal equilibrium. However, as shown in [l5[, a dynamical 
violation of CPT (which implies a different spectrum of particles and antiparticles) may give 
rise to the baryon number asymmetry also in a regime of thermal equilibrium. A successful 
mechanism for explaining the asymmetry between matter and anti-matter is provided by 
Leptogenesis [l^. In this scenario, in which the Majorana neutrino is introduced, it is 
possible to generate the baryon asymmetry if the asymmetry is generated in the lepton 
sector at either GUT or intermediate scales, even if the baryon number is conserved at high 
energy scales. 

In the present paper we show that the coupling of the heavy neutrinos with the gravita- 
tional background gives rise to an effective potential that modifies differently the dispersion 
relations of neutrinos with left and right helicity, leading to a net lepton asymmetry even if 
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neutrino are in thermal equilibrium^. The lepton asymmetry is then converted into baryon 
asymmetry by the action of sphalerons in the elctroweak era. For this mechanism to work, a 
time varying non-zero Ricci curvature is required during the radiation dominated era. This 
is provided by f{R) cosmology (mechanisms for baryo/leptogenesis, similar to one discussed 
in this paper, can be found in [17|). 

The paper is organized as follows. In Section II we discuss the leptogenesis mechanism 
based on the curvature coupling of heavy neutrinos. In Section III we derive the main 
equations of f{R) models of gravity, studying also the constraints provided by BBN physics. 
In Section IV and V we investigate the baryon asymmetry for some specific models of f{R)- 
Conclusions are shortly discussed in Section VI. 



2. CP ODD GRAVITATIONAL COUPLING OF MAJORANA NEUTRINOS 



In this Section, we study the generalization in the matter Lagrangian by including higher 
order terms in R consistent with general covariance, Lorentz-invariance in a locally iner- 
tial frame. Consider the action for a four component Dirac fermion ip which couples to 
background gravity, 

S^[g,.,^] = j d'xV^[t^Y(^,-^,)^- (2.3) 

-/ll(i?)V^V^-2/l2(i?)^75^] , 

where hi{R) and h2{R) real valued scalar functions of the curvature, 

h{R) = M + g,{R), 

h^iR) = M' + g^iR). (2.4) 

Here hi is a generalization of the neutrino mass term. Note that since ipj^ip transforms as 
a pseudo-scalar, the h2 term is odd under CP. We write the four- component fermion 

The lagrangian in terms of the two-component fields ipR and ipi becomes 

- /ii(^i^L + ^i^i?) - I /i2(V'i^L - ^i^i?) , (2.6) 

where = (/, a*) and = (/, — o"*) in terms of the Pauli matrices. The /12 term can be 
rotated away by a chiral transformation 

^i^e-*"/VL ^ii^e^°/Vi?- (2.7) 



^ The coupling of the chiral fields to the cosmological background is odd under CP which biases the gen- 
eration of leptons over anti-leptons in the presence of lepton number violating interactions at thermal 
equilibrium (similar to [l3|). The thermodynamic interpretation of such a scenario is described, as usual 



15, 



ITj , by the introduction of an effective chemical potential which is different for particles and antipar- 



ticles owing to the CP odd coupling. 
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Keeping terms to the linear order in a, we see that the lagrangian (12. 6 p changes by the 
amount 

- hi{ia) (i)[i)R - ^i^i) - ih2{ia) (^^{^R + V'LV'l) • (2.8) 

Now we choose a = — /i2/^i to ehminate the chiral mass term and obtain for the total 
Lagrangian 



C = ^^^^a>'(^^-'^^)iJR + li;ia^'(^^-'^^)^Pl 



- ;^(/^? + /^2)(4^h + <^^l) . (2.9) 

If hi and h2 are constants then, one can always rotate the axial-mass term away. We will 
assume that the neutrino mass M ^ gi therefore /ii ~ M and since a constant M' can 
be rotated away /12 = 92- Further we will assume that the background curvature is only 
dependent on time. The lagrangian ( 12. 9p then reduces to the form 

C = a^(^,, - t,)tljR + a^i^^ - ^,)tPL - 

_,k(|),,.,„(||)_ 

- M(^^[^R + ^ltPL) . (2.10) 
The equation of motion for the left and the right helicity fermions derived from (I2.10p are 

ta^d^i^L - (^] - MijR = . (2.11) 



Written in momentum space ip{x) = ■?/)(p)e*(-^*~^'^) the equation of motion of ipR and ipL are 



EL + ^^iJL + a-pijL-MtPR = 0. (2.12) 
The canonical momenta of the ipL and ipR fields are as usual 

TTl = = t^l, TTR = = zV^jj . (2.13) 
Then the canonical Hamiltonian density is 

= iijj[ipL + iip[(T ■ VipL + ii^R^JR - iipRO- ■ VijjR + M (ij[tljR + tpl^ipL 



where we have introduced the number density operators of the left and right chirahty modes, 



riL = V^i^L, nR = ip]^ipR . (2.15) 
The partition function in terms of this effective Hamiltonian is 

Z = Tre'^^ = Tre"'^(^°~^^"^""«''«) , (2.16) 

where /3 = 1/T and Ho is the free particle Hamiltonian. We see that when g2 is non-zero 
then the effective chemical potential for the left chirality neutrinos is hl = —g2/M and for 
the right-chirality neutrinos is fiR = 5'2/M. In the presence of interactions which change 
ipL "ipR at thermal equilibrium there will be a net difference between the left and the right 
chirality particles, 
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j d^p 
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M ■ 
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X _|_ qP{Ep-hr) I _|_ qI3{Ep-iil) 



(2.17) 



In this paper we consider the simplest case of h2{R) = g2{R) = R/mp as a linear function 
of the curvature R. The axial term in f l2.3p is a CP violating interaction between fermions 



and the Ricci curvature described by the dimension-five operator [18|, 119 



mp 

This operator is invariant under Local Lorentz transformation and is even under C and odd 
under P and conserves CPT. In a non-zero background R, there is an effective CPT viola- 
tion for the fermions. Take ip = {Np, Np)'^, where Np is a heavy right handed neutrino and 
Np a left handed heavy neutrino, which decay into the light neutrinos. Majorana neutrino 
interactions with the light neutrinos and Higgs relevant for leptogenesis, are described by 
the lagrangian 

£ = -/i„^(0t W^Il/,) - ^N'pMNp + h.c. , (2.19) 

where M is the right handed neutrino mass-matrix, lia = {^a,€.~)J^ is the left-handed 
lepton doublet (a denotes the generation), = (0"'",0°)"^ is the Higgs doublet. In the 
scenario of leptogenesis introduced by Fukugita and Yanagida, lepton number violation is 
achieved by the decays Np cj) + II and also Np^ — 0^ + l^^. The difference in the 
production rate of II compared to Z£, which is necessary for leptogenesis, is achieved via 
the CP violation. In the standard scenario, n{Np) = nlNf^) as demanded by CPT, but 
T{Np — > /l + 0) 7^ ^{^R ^ + 0^) due to the complex phases of the Yukawa coupling 
matrix /iq,^, and a net lepton number arises from the interference terms of the tree-level and 
one loop diagrams j20|, 2ll . 



In our leptogenesis scenario we have T{Np — )■ /l + 0) = ^{^p + 0^) but there is a 
difference between the heavy light and left chirality neutrinos at thermal equilibrium due to 
the CP violating gravitational interaction fl2.18p . 

n(iVH)-n(iVa = y^. (2.20) 
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The Nr -h- interaction can be achieved by the scattering with a Higgs field. For example 
in SO(IO) theories Majorana masses for the right-handed neutrinos are generated either 
from 

16f 16f T26h DyS'N'^Nn, (2.21) 
or from the non-renormalizable operators suppressed by some mass scale A 

^16fl6fT6HT6H3 {^'iV^iV^. (2.22) 

When the GUT Higgs fields S' or S acquire a vev, a large Majorana mass M is generated 
for Np which breaks lepton number spontaneously. This following the see-saw mechanism 
leads to small neutrino masses at low energies. At temperatures larger than the heavy 
neutrinos and the GUT Higgs masses one there will be chirality fiip scattering interactions 
like S + Nji i-^ S + N"^ which change the lepton number. The interaction rate is given by 

T{SNr<^ SN'^) = M = ^ (J^^ T\ (2.23) 

The interactions decouple at a temperature when r(TD) = H{To) from which we derive 
the decoupling temperature to be 

To = 13.7. (^)' ^ = {^y^^, (2.24) 

where we have used M = f{S)/A. The lepton number asymmetry is frozen in the heavy 
neutrino number asymmetry at temperature To and is passed on to the light lepton sector 
after the subsequent decays of into light neutrino. Substituting (12.24^ in the expression 
for lepton asymmetry fl2.20p we obtain the value of frozen in lepton asymmetry as 

.„ ^ "<^"' - = 4 (2.25) 

where R is to be evaluated at T = Tq. In the subsequent sections we first give a survey of 
f{R) gravity and then compute R in f{R) cosmology. 

A comment is in order. In the case in which the fermion is, for example, an electron 
one also gets a splitting of energy levels E^cr) — E{eL), but this does not lead to lepton 
generation of lepton asymmetry as both ei and ep carry the same lepton number. 



3. FUNDAMENTAL EQUATIONS AND BBN IN f{R) THEORIES OF GRAVITY 

As pointed out by Nojiri and Odintsov in [5], modified gravity contains many topics which 
make these models very attractive, as for example: 1) they provide a natural unification of 
the early-time infiation and the later-time acceleration of the Universe owing to the different 
role of the gravitational terms relevant at small and large scales; 2) they allow to unify dark 
matter and dark energy; 3) they provide a framework for the explanation of hierarchy 
problem and unification of GUT with gravity. However, many f{R) models of gravity are 
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strongly constrained (or ruled out) by solar system tests. In this respect, available models 
are those proposed in Ref. [2^ 



f{R) = -m 



2 ci(i?Vm2)2- 



l + C2{R/m 



2\2n ' 



and in Ref. 23 



f{R) =R + XR. 



St 



R^ 

Rlt 



- 1 



(3.26) 



(3.27) 



where A, m^, Ci^2; Rst-, n, and d are free parameters. 

These functions can be expanded in the appropriate regimes, reproducing simplest form 
of f{R). A particular subclass is given by 



f{R) =R + aR'' 



(3.28) 



where a > has the dimensions [energy] ~2("~^) and n > 0. Combinations of the free 
parameters allow to get a description of cosmic acceleration. Particularly interesting is the 
case n = 2 (referred in literature as Starobinsky's model [23]) 



f{R) =R + aR\ 



(3.29) 



The model fl3.29p has been studied in the framework of astrophysics and cosmology. For 
instance, gravitational radiation emitted by isolated system constraints the free parameter 
to |a| < (10^^ — 10^*)m2 25|, |26|. Eot-Wash experiments lead instead to the constraints 



lal < 2 X IQ-^m^ . 



(3.30) 



However, more stringent constraints are provided by the Cosmic Microwave Background 
(CMB) physics. The amplitude of the curvature perturbation corresponding to f l3.29p is 



Pi 



~ — ^ ^, with Nk ~ 55. Using the WMAP 5-years data (Pn ~ 2.445 x lO'^) HJ, it 

ISvr arrip 



follows that a is constrained as 



lal < lO^^W 



(3.31) 



The bound f l3.3ip is obtained in the regime R ^ (in this regime the model describes 
the inflationary epoch). 

Finally, we mention the model 

f{R)=R'+\ (3.32) 

where e ^ 1. The tightest bound on e has been inferred by Clifton and Barrow 28| in the 
framework of perihelion precession of Mercury 



< e < 7.2 X 10 



-19 



(3.33) 



The field equation obtained by the variation of the action (11. ip with respect to the metric 



are 



(3.34) 
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where the prime stands for the derivative with respect to R. The trace reads 

3Df' + f'R-2f = KX, (3-35) 

with Tg the trace of the energy-momentum tensor. 

In the spatially flat Friedman-Robertson- Walker (FRW) metric 

ds'^ = -dt^ + a^{t) [dx^ + + dz^] , (3.36) 

Eqs. fl3:Ml) and fl335|) become 

3f'H^ - ^ + 3Hf"R = n^p , (3.37) 

-2fH^ - fie + nHf'R -R) = k\p + p), (3.38) 
3/'"i?2 + 3f"R + QHf'R + f'R - 2/ = k^t^ ^ (3.39) 

where H = d/a, T°q = — p, T^^- = and Tg = p — 3p. Moreover, the Bianchi identities 
give a further condition on the conservation of the energy 

p + 3-{p + p) = 0. (3.40) 

a 

In what follows, we shall look for those solutions of field equations such that the scale factor 
evolves as 

a{t) = aot^ , H = -. (3.41) 

t 

The scalar curvature turns out to be 

R = 6{2H' + H) = . (3.42) 

The f{R) model we concern here is that one of Eq. (13.281) . By using Eqs. fl3.37p and fl3.38p 
and the usual expression relating the energy density and the pressure, p = wp, where w is 
the adiabatic index, one gets 

w = ^ + S{t) , (3.43) 

where 

2 f(3 + nA 



with 

A = aR''-^[(3{2 -n)-{n- l){2n - 1)] . 
The energy density p assumes the form 

«V^f(l + |). (3.45) 

Notice that during the radiation dominated era (/3 = 1/2), to which we are mainly interested, 
the quantity A vanishes because i? = 0, as well as the perturbation 6, and the adiabatic 
index reduces to the standard value w = 1/3. 
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For our estimations on the baryon asymmetry, it is relevant to consider the following two 
features. 

Let us first discuss the BBN constraints on f{R) theories of gravity (see also [2^ 



29|). 



For the BBN scenario, one has to consider the weak interaction rate of particles {p, n, 
and v) in thermal equilibrium. For T ^ Q {Q = m. 



nip, where Tn„p are the neutron 



and proton masses), one gets A(T) ~ gT^, where g = 9.6 x 10 "^^eV ^ 0, 30[. To estimate 



where A 



e-(«--*/)/^. tf 



the primordial mass fraction of '^He, one defines Yp = X i^^f^^^-^ 
and tn are the time of the freeze-out of the weak interactions and of the nucleosynthesis, 
respectively, r ~ 887sec is the neutron mean life, and x{tf) = e~^^'^^^f^ is the neutron to 
proton equilibrium ratio. The function A(t/) represents the fraction of neutrons that decay 
into protons in the time t G Deviations from Yp (generated by the variation of the 

freezing temperature Tf) are given by [sij 



SY, = Y„ 



Y„ 



2A 



2A 



^ In ;-: - 1 - z:^ 



Y„ 



2t 



T 



ST 



f 



f 



(3.46) 



In the above equation we have set 5T{tn) 
energy. The current estimation on 1^, Yp 

6Tf 



because Tn is fixed by the deuterium binding 



0.2476 ± 6Yp, with \6Yp\ < 10"^ |32|, leads to 



T 



f 



< 4.7 X 10^ 



(3.47) 



The freeze-out temperature T is determined by A 
Tf ~ 0.6 MeV and 



H. One gets T = Tf{l + f^), where 



ST 



f 



6 



Tf 15 
Equations IK^ and (K^ imply 



5 qm pTf 



1.0024 (3 - 



2/3 - 1 < 9.4 X 10" 



(3.48) 



(3.49) 



As we shall see, the observed value of ^^exp, Eq- (ll-2p . is obtained in our model of Leptoge- 
nesis for values of 2/3 — 1 well below the bound f l3.49p . 

The second aspect we wish to discuss concern the regime for which the terms that modify 
general relativity induce very tiny deviations from the standard cosmological model, i.e. 



a 



BJ"'^ < 1 . 



The latter implies that the parameter a is bound from above by 

a 



m 



2\n-l 



< n, 



where 



n = 



1 



45 



[6/3|2/3- l|]"-i \16n^g. 



n-l 



-^Q_68(n-1) .^^x 4(n-l) 



In this approximation {A <^ /3) one also obtains 

(n - 1)^ 



. 2 
~ — 

3/3 



2/3 



/3 



< 1. 



(3.50) 



(3.51) 



(3.52) 



(3.53) 
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4. GRAVITATIONAL LEPTOGENESIS AND f{R) GRAVITY 



As already noted, in the standard cosmological model R vanishes during the radiation 
era. However, tiny deviations from General Relativity give rise to a Ricci curvature, as 
well as its first time derivative, different from zero so that a net lepton asymmetry can be 
generated. By making use of the definition of Ricci scalar curvature f l3.42p . it follows 

R^JimzA, (4.54) 

Using the relation between the cosmic time and the temperature - ^ ( ^^^/^ y. 

t 90 mp 

(9 , the parameter characterizing the lepton asymmetry (12.25^ assumes the form (to 
leading order in (2/3 — 1)) 

.(2/; -1)3.4X10--^°"°^^^ 



M VlOiSGeV 

In this equation To is the decoupling temperature at which the lepton number violating 
GUT reactions which change the number of and N'^ go out of equilibrium. This lepton 
asymmetry can be converted into baryon asymmetry by the action of sphalerons in the 



electroweak era 01. 



An inspection of f l4.55p immediately revels that the observed baryon asymmetry can be 
obtained for ~ lO^^^GeV, M ~ lO^GeV, provided that 2/3 - 1 ~ 3 x 10"^ 

Another possibility to get rj ~ 10^^*^ is given by taking, for example, {Tp, ~ lO^^GeV; 
M = lO^^QgV, 2/3 - 1 ~ 10-15 - 10-13} or {Tp ~ lO^^GeV; M = lO^GeV, 2/3 - 1 ~ 
XQ-is _ I0~i^}. The value of the heavy neutrino mass ~ lO^^GeV is consistent with the 
atmospheric neutrino scale = 0.05 eV, obtained from the see-saw relation rrii, = mj^/M 
with the Dirac mass scale mp ~ (9(10) GeV. Similar conclusions holds also for M = lO^GeV. 

The lepton asymmetry generated via fl4.55p is passed on to the light neutrino sector 

when the heavy neutrino decays at temperature T ~ M ~ (10^ — 10i^)GeV. The effects of 

washed out are avoided by considering the effective (five dimensional) operator violating the 

(J 

lepton number AL = 2, Cw = qj^^^^^ 0*)(0^ lp) + h.c.. We suppress the generation indices. 

Ip = {u, e~)\ is the left-handed lepton doublet, (p is the Higgs doublet, and M corresponds in 
general to some large mass scale (identified in our case with the heavy neutrino mass) j33 |. 

122 ICPT^ 

The interaction rate is T{i'p + </> o z//j + 0) = {n^ a) = electroweak 

era, when the Higgs field in Cw acquires a vev, {(p) = v = 174GeV, the five dimensional 
Weinberg operator gives rise to a neutrino mass matrix = ^jj-. The lepton number 
violating interactions decouple when T{Ti) = H{Ti), which implies 

r, = (2ffliM^^2xlO"f^VGeV. (4.56) 

ml nip \ my J 
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TABLE I: Estimations of 11 for varying n. The values of and 2/3 — 1 are those leading to 
rj ~ 10-1°. We fix Td = lO^^GeV (similar results follow for Td = lO^'^GeV). On the left are 
reported the values for 2/3 — 1 = 10^^^ on the right for 2/3 — 1 = 10"^^. 



n 



n 







1.1 


10 


-6 


1.1 10-6 


1.5 


10" 


-27 


1.5 10-24 


3 


10" 


104 


3 10-92 



n 



n 



Therefore the heavy neutrino decays occur at a temperature T c:^ M c:^ (lO^ — 10^2) GeV well 
below the temperature T; = 2 x lO^^GeV at which the light-neutrino lepton number violating 
interactions are effective. As a consequence, the lepton number asymmetry from the decay 
of asymmetric number of heavy neutrino decays is not washed out by Higgs scattering with 
light neutrinos. 



4.1. The general case n ^ 

Here we discuss the model (I3.28P with generic n. In Table I is reported the estimation of 
n given by fl3.52p for different values of n, Tjj and 2/3 — 1. We see that 11 <^ 1 for n > 1. 
The values of the decoupling temperature and 2/3 — 1 are those before used for inferring 
the correct estimation of rj. By making use of fl3.50p one can properly fix a in order that 
our approximation holds. 

The case n < 1 deserves also a discussion. As shown in jsf, this case yields the constraint 
n < IQ-io if considered as a available candidate for dark matter (the authors assume a < 0). 
Such a bound is obtained from the combination of solar system experiments with the the 
equivalence principle violation experiments. Using f l3.52p with = lO^^GeV and 2/3 — 1 ~ 
10^22^ Qj^g obtains IT ~ 10^^ order that our approximation work, it is requires a parameter 
a bounded at least by a < 10~^^{m'^)"-~^ . Also in this model a net lepton asymmetry of the 
order of ~ 10~^^ can be therefore obtained. 

We do not investigate the case n < because these f{R) models of gravity are affected 
by instability problems 



4.2. The case n = 2 



A particular subclass of the model studied in the previous Section is the case in which 
n = 2. The f{R) function is given by (13.290 . The parameter a now has the dimensions 
m2 ([energy] -2)^ As already pointed out in the Introduction, Eot-Wash experiments provide 
stringent constraints on a given by (I3.30p . On the other hand, more stringent bound can 
be inferred from CMB physics (I3.3ip . 

For n = 2, the function 11 reads 

n = 1.8x10-61— . 

2/3 - 1 V T J 
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Taking 2/3 - 1 = lO'^ and To = lO^^GeV one gets H ~ 10"^^ while for 2/3 - 1 = 10"^^ and 
Td = lO^^GeV it follows 11 ~ 10"'^^. These values are well below the upper bounds fl3.30p 
and (I33T]) . In both before discussed, the observed baryon asymmetry is obtained. 



5. THE MODEL f{R) = R^+^ 

In this Section we finally study the model 




(5.57) 



where A is a constant with dimensions m^p^^^^\ Eqs fl3.4ip . (13.371) and (I3.38P imply 

l + e = 2/3. (5.58) 

The constraint on the parameter e given by (13.331) and Eq. (15.581) yield [2^ 

< 2/3 - 1 < 7.2 X 10"^^ (5.59) 

Taking 2/3 - 1 ~ 7 x 10^^^ M ~ lO^GeV Q and Td ~ 5 x lO^^GeV, the lepton asymmetry 
given by Eq. (I4.55P yields rji ~ 10~^°. These results show that this model of extended 
theory of gravity can work for our mechanism. 



6. CONCLUSIONS 

In this paper we have shown that f{R) models of gravity, described by (I3.28p . (13.290 . 
and (13.320 . are available frameworks for the generation of a net baryon asymmetry. These 
models, besides to be able to give a dynamical description of dark energy and to be also 
favorite candidate for the dark matter puzzle, are consistent with solar system experiments. 

The matter-antimatter asymmetry scenario discussed in this paper relies on a CP violat- 
ing gravitational interaction between the heavy Majorana neutrinos and the Ricci curvature. 
The lepton asymmetry generated in this way is in thermal equilibrium during the radiation 
era at GUT scale. The subsequent decays of these heavy neutrinos into the light stan- 
dard model particles and the conversion (via sphaleron processes) of lepton asymmetry into 
baryon asymmetry can explain the observed baryon asymmetry of the universe. As we have 
shown, the leptogenesis mechanism works for decoupling temperatures of the order of GUT 
scale, Td ~ (10^^ - 10^^)GeV, heavy neutrino mass M ~ (10^ - W^)GeY, and finally for 
very tiny deviations from the standard evolution of the Universe during the radiation era. 
I.e. a ~ with /3 < 1/2 + 10"^ (a bound below the upper bound obtained from BBN 
2/3 - 1< 10-3). 
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